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ABSTRACT 


Free  Electron  Lasers  are  the  focus  of  a  recently  announce  Innovative  Navy 
Prototype  to  develop  a  directed  energy  weapon  system  for  the  self-defense  of  ships. 
Operating  in  a  shipboard  environment  poses  several  challenges  that  must  be  overcome. 
Short  Rayleigh  length  systems  offer  solutions  to  some  of  these  problems.  Simulations 
were  performed  to  examine  the  benefit  of  short  Rayleigh  length  designs  in  the  face  of 
electron  beam  misalignment.  Additionally,  simulations  were  performed  to  explore  the 
effect  of  quadrupole  misalignment  on  electron  beam  position  and  trajectory,  and 
ultimately  on  FEL  performance. 
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I.  INTRODUCTION 


The  United  States  Navy  has  a  long  history  of  embracing  new  technologies  for  the 
purpose  of  furthering  control  of  the  seas  and  littoral  areas  of  the  world.  The  Navy  is 
continuing  this  tradition  by  funding  research  for  directed  energy  weapons.  High-energy 
laser  systems,  which  have  the  potential  to  fundamentally  change  the  way  surface  vessels 
operate,  are  an  important  area  of  this  research.  A  high-energy  laser  could  engage  a  wide 
variety  of  targets,  including  anti-ship  missiles,  theater  ballistic  missiles,  aircraft,  other 
surface  vessels,  and  potentially  even  near-surface  mines.  Equation  Chapter  1  Section  1 

High-energy  lasers  have  already  been  developed  and  tested  for  military  purposes. 
A  joint  US  Army  and  Israeli  project  has  resulted  in  the  development  of  the  tactical  high- 
energy  laser  (THEL),  a  mid-infrared  advanced  chemical  laser  (MIRACL)  that  has 
successfully  destroyed  Katyusha  rockets  in  flight  [1],  The  US  Air  Force  has  developed  a 
chemical-oxygen-iodine  laser  (COIL)  for  use  in  Air  Borne  Laser  (ABL)  program  [1]. 
The  US  Navy  recently  approved  an  Innovative  Navy  Prototype  (INP)  to  build  a  Free 
Electron  Laser  (FEL)  with  the  goal  of  reaching  the  100  kW  power  level  milestone. 

Unlike  chemical  lasers,  an  FEL  makes  use  of  a  relativistic  beam  of  electrons 
passing  through  a  periodic  magnetic  field  to  produce  coherent  electromagnetic  radiation. 
The  components  necessary  for  this  method  of  lasing  make  FELs  larger,  heavier,  and  more 
expensive  than  many  other  lasers  used  for  low  power  applications.  However,  the  FEL’s 
high  power  availability  and  tunable  wavelength  make  it  very  useful  in  a  wide  range  of 
applications  at  both  low  and  high  power.  FELs  also  scale  up  very  favorably,  meaning 
that  significantly  more  powerful  FELs  need  not  be  very  much  larger  or  more  costly  than 
lower  power  versions.  Additionally,  an  FEL  is  an  all-electric  system  that  could  run  off  of 
ship’s  power  without  the  need  for  refueling  of  dangerous  chemicals  or  explosives,  or 
venting  of  toxic  exhaust  fumes.  Figure  1  is  an  illustration  of  how  the  Navy  might  employ 
an  FEL  as  a  weapon  system. 
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Figure  1 .  Concept  drawing  of  a  shipboard  FEL  weapon  system. 


This  thesis  investigates  the  effects  of  quadrupole  magnet  misalignment  on  FEL 
operational  performance.  This  study  builds  upon  work  done  previously  ([2])  to 
investigate  the  effects  of  shifting  a  single  quadrupole,  and  is  aimed  at  establishing 
vibration  dampening  requirements  given  the  effect  of  shifting  multiple  quadrupoles. 
Chapter  II  provides  some  background  information  about  the  FEL’s  development  and  its 
suitability  as  a  naval  weapon  system.  Chapter  III  describes  the  major  components  of  an 
FEL.  Chapter  IV  discusses  the  physics  of  FEL  operation.  Chapter  V  describes  how  the 
magnetic  elements  of  the  electron  beam  transport  system  affect  electron  trajectories. 
Chapter  VI  describes  simulations  that  were  conducted  to  detennine  the  affects  of 
quadrupole  magnet  misalignment.  Finally,  Chapter  VII  provides  an  analysis  and 
summary  of  the  simulations,  and  suggests  directions  of  further  research. 
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II.  BACKGROUND 


It  is  a  basic  principal  of  electrodynamics  that  a  charge  undergoing  an  acceleration 
will  radiate  electromagnetic  waves  [3],  This  principal  often  factors  into  the  operation  of 
particle  accelerators,  which  use  electric  and  magnetic  fields  to  accelerate  charged 
particles.  The  development  of  particle  accelerators,  as  well  as  numerous  and  varied  other 
devices  that  make  use  of  this  principal,  is  an  important  part  of  the  history  of  the  FEL. 
What  follows  is  not  an  exhaustive  history,  but  rather  a  broad  overview  of  significant 
milestones  in  the  development  of  FELs.  First,  however,  it  is  instructive  to  discuss  a  few 
related  devices.  Equation  Chapter  2  Section  1 

A.  RELATED  DEVICES 

1.  Magnetron 

Magnetrons  are  vacuum  tubes  in  which  magnetic  fields  are  used  to  cause 
electrons  from  a  cathode  to  spiral  around  resonator  cavities  in  order  to  generate  coherent 
microwave  radiation.  The  first  magnetron  was  developed  by  Albert  W.  Hull  of  General 
Electric  in  the  early  1920’s.  Magnetrons  have  been  used  as  amplifiers  in  radio  receivers, 
and  as  sources  in  various  radar  systems.  They  are  capable  of  operating  at  very  high 
efficiency,  around  60%.  Today  they  can  be  found  in  a  wide  range  of  applications, 
including  microwave  ovens  in  the  home  [4]. 

2.  Cyclotron 

Invented  in  1929  by  Ernest  Lawrence  of  UC  Berkley,  the  cyclotron  is  a  particle 
accelerator  that  offered  significant  advantages  over  linear  accelerators  (linacs)  previously 
in  use.  Accelerators  work  by  using  radio  frequency  (RF)  waves  to  impart  energy  to 
charged  particles.  The  size  and  cost  of  linacs  placed  an  upper  limit  on  the  energies  that 
could  practically  be  achieved.  Cyclotrons  use  a  magnetic  field  to  circulate  charged 
particles,  making  them  much  smaller  and  much  more  efficient  then  the  linear  equivalent. 
The  upper  limit  on  energy  levels  achievable  with  a  cyclotron  results  from  relativistic 

3 


effects  that  cause  the  particle’s  mass  to  increase,  which  disrupts  the  particle’s  orbit  and 
causes  it  to  move  out  of  phase  with  the  RF  waves  [5], 

3.  Klystron 

The  klystron,  invented  by  the  brothers  Russell  and  Sigurd  Varian  of  Stanford 
University  in  1937,  is  a  device  that  uses  an  evacuated  electron  tube  and  two  resonant 
cavities  to  amplify  a  reference  signal  of  microwave  or  radio  frequency  radiation  [6]. 
Klystrons  can  be  used  to  produce  low  power  reference  signals  for  radar  systems,  higher 
power  carrier  waves  for  communications,  and  very  high  power  driving  waves  for  particle 
accelerators  They  can  range  in  frequency  from  100’s  of  MHz  to  100’s  of  GHz,  and 
produce  power  levels  up  to  100’s  of  MW. 

4.  Synchrotron 

The  synchrotron  was  independently  invented  by  Edwin  McMillan  of  UC  Berkley 
and  Vladimir  I.  Veksler  of  the  Soviet  Union  around  the  same  time  in  1945  [5].  It  was 
conceived  as  a  particle  accelerator  that  uses  a  magnetic  field  to  circulate  charged 
particles,  and  an  electric  field  to  accelerate  the  particles.  The  device  takes  its  name  from 
the  fact  that  the  electric  and  magnetic  fields  must  be  synchronized  with  the  particle  beam. 
This  synchronization  takes  advantage  of  a  so-called  “phase  stability”  to  avoid  the 
problems  of  a  cyclotron  at  relativistic  energies  [7].  In  addition  to  accelerating  electrons, 
synchrotrons  produce  radiation.  At  first  this  side  effect  was  considered  a  nuisance,  since 
it  represented  an  unwanted  loss  of  energy.  It  has  since  been  recognized  that  synchrotron 
radiation  possesses  some  interesting  characteristics,  such  as  high  polarization,  high 
brilliance,  and  low  emittance.  These  characteristics  make  synchrotron  radiation  useful 
for  a  wide  range  of  research  and  industrial  applications  [8].  The  tenn  synchrotron 
radiation  now  applies  to  all  radiation  from  accelerating  charged  particles  of  relativistic 
energy,  whether  this  takes  place  in  a  synchrotron  or  in  some  other  device. 
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B.  HISTORY  OF  FEL  DEVELOPMENT 

Since  FELs  make  use  of  synchrotron  radiation,  the  first  observation  of  this 
radiation  is  a  logical  place  to  begin  the  history  of  FEL  development.  This  observation 
was  made  in  1947  by  researchers  using  a  General  Electric  70  MeV  synchrotron 
accelerator  [9],  The  discovery  of  this  radiation  was  not  the  goal  of  the  General  Electric 
researchers,  nor  was  it  a  surprise.  The  presence  of  this  radiation  had  already  been 
predicted,  and  it  was  speculated  that  it  would  place  an  upper  limit  on  the  energy  levels 
attainable  by  circular  accelerators  [10]. 

One  of  the  earliest  ancestors  of  the  FEL  is  described  in  a  1952  patent  by  Elmer  J. 
Gom  [11],  [12],  The  patent,  which  was  submitted  in  1947,  is  for  a  “Traveling-Wave 
Electron  Reaction  Device.”  The  stated  purpose  of  this  device  was  to  cause  electrons  to 
interact  with  an  electromagnetic  wave  in  such  a  way  as  to  either  give  up  energy  to  or 
extract  energy  from  the  wave.  In  the  first  case,  the  device  would  be  an  amplifier,  and  in 
the  second  case  it  would  be  an  electron  accelerator.  This  device  was  certainly  no  FEL,  as 
the  electron  beam  was  not  relativistic,  and  the  emitted  light  was  not  very  coherent. 
However,  it  did  employ  the  concept  of  an  electron-optical  wave  interaction  that  relied  on 
a  periodicity  in  the  electron  beam. 

In  the  1950’s,  H.  Motz  of  Stanford  University  explored  the  radiation  from  fast 
electron  beams  passing  through  a  series  of  magnetic  fields  with  alternating  polarity, 
called  an  undulator  [13-16].  It  was  discovered  that  undulators  could  produce 
synchrotron  radiation  more  efficiently  than  synchrotron  bending  magnets  could.  Using  a 
3  MeV  electron  accelerator,  Motz  produced  millimeter  wave  radiation.  Later,  using  a 
100  MeV  accelerator  he  was  able  to  produce  visible  light.  The  radiation  he  produced  was 
not  very  coherent,  though  Motz  did  discuss  the  possibility  that  more  coherent  light  might 
be  produced  by  this  method.  His  work  with  undulators  proved  the  versatility  of  this 
device  by  demonstrating  operation  over  a  range  of  wavelengths.  Undulators  are  now 
used  in  conjunction  with  storage  rings  to  produce  synchrotron  radiation  for  a  wide  range 
of  research  and  industrial  applications  [5]. 
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Another  important  development,  made  by  R.  M.  Phillips  in  the  late  1950’s,  was 
the  ubitron.  Somewhat  similar  to  a  klystron,  this  device  takes  its  name  from  the  phrase 
“undulated  beam  interaction”  [17].  The  ubitron  produced  centimeter  and  millimeter 
wavelength  radiation  at  much  higher  powers  than  klystrons  or  magnetrons  of  the  time, 
and  was  once  evaluated  for  its  potential  as  a  directed  energy  weapon  [11]. 

The  devices  described  so  far  either  used  resonator  cavities  to  produce  coherent 
microwave  radiation,  or  had  no  resonator  mechanism  and  produced  incoherent 
synchrotron  radiation.  FELs  use  either  an  open  resonator  or  a  seed  laser  to  create 
coherent  emission  of  optical  range  synchrotron  radiation.  The  first  person  to  demonstrate 
the  operation  above  threshold  of  a  true  FEL  was  J.  H.  J.  Madey,  of  Stanford  University 
[18].  In  1977,  a  free  electron  laser  in  oscillator  configuration  was  operated  above 
threshold  at  a  wavelength  of  3.4  microns. 

C.  APPLICATIONS  OF  A  FREE  ELECTRON  LASER 

FELs  have  a  great  deal  of  promise  in  a  wide  range  of  applications.  One  of  the  key 
advantages  of  an  FEL  is  wavelength  “tunability”,  which  makes  this  type  of  laser 
extremely  versatile.  As  will  be  explained  in  more  detail  later,  the  output  wavelength  of 
the  laser  can  be  selected  in  the  design  process,  and  can  be  adjusted  even  after  the  system 
has  been  constructed. 

This  “tunability”  makes  the  FEL  a  good  candidate  for  improving  manufacturing 
in  the  areas  of  polymer  surface  processing,  micromachining,  and  metal  surface 
processing,  to  name  a  few  [19].  The  university  of  Hawaii  is  developing  an  FEL  as  a  light 
detection  and  ranging  (lidar)  transmitter  for  remote  sensing  applications  [20]. 
Additionally,  FELs  offer  many  advantages  in  the  medical  field.  Their  selectable 
wavelength  would  make  possible  a  reduction  in  scattering,  which  would  improve 
resolution  in  imaging  applications.  It  may  also  be  possible  to  reduce  collateral  damage  in 
cutting  applications  by  reducing  absorption  and  heating  in  un-targeted  tissues  [21]. 
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D. 


FEL  AS  A  WEAPON 


A  laser  causes  damage  to  a  target  by  delivering  a  large  amount  of  energy  to  a 
small  area,  causing  a  portion  of  the  target  to  reach  its  melting  temperature.  Consider  a 
sample  of  aluminum,  10  cm  by  10  cm,  and  1  cm  thick.  Aluminum  is  a  common  material 
in  the  construction  of  missiles,  and  the  destruction  of  this  size  sample  should  be  enough 
to  cause  a  missile  to  break  apart  in  flight  [1].  The  physical  properties  of  aluminum  are 
summarized  in  Table  1. 


Density 

2700  kg/nfi 

Thermal  conductivity 

237  W/m  K 

Specific  heat  capacity 

900  J/kg  K 

Melting  temperature 

933.47  K 

Emissivity 

0.5 

Heat  of  fusion 

1 1  kJ/mol 

Table  1.  Physical  properties  of  Aluminum,  from  [22] 


The  power  radiated  away  by  our  heated  sample  of  aluminum  is  given  by 

P  =  eaA(T*-T^)t  (2.1) 

where  s  is  the  emissivity,  a,  is  the  Stefan-Boltzman  constant,  A  is  the  area,  TAI  is  the 
temperature  of  the  aluminum  and  Tenv  is  the  temperature  of  the  surrounding  environment. 
The  power  conducted  away  is  given  by 

H  =  kA(TAl  -  Tenv)  /  Ax ,  (2.2) 

where  k  is  the  thermal  conductivity,  and  At  is  the  thickness.  Our  sample  has  a  surface 
area  of  100  cm  ,  and  a  thickness  of  1  cm.  Taking  Tenv  to  be  300  K  and  TAl  to  be  the 
melting  point  of  aluminum,  (2.1)  yields  P  =  214  W  radiated  away,  while  (2.2)  yields  H 
=  150  kW  conducted  through  the  back  of  the  sample  and  H  =  60  kW  conducted  through 
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the  edges.  The  total  power  leaving  our  sample  of  aluminum  when  it  is  near  its  melting 
point  is  a  little  over  200  kW.  If  more  than  this  amount  of  power  is  applied  to  the  sample, 
the  aluminum  will  melt. 

From  the  density,  we  find  our  sample  is  about  0.3  kg.  The  heat  energy  required  to 
bring  this  amount  of  aluminum  from  300  K  up  to  the  melting  temperature  is  given  by 

Q  =  cQm(TAI  -  Tem) ,  (2.3) 

where  c0  is  the  specific  heat  capacity  and  m  is  the  mass.  Substituting  from  Table  1 
results  in  Q  =175  kJ.  One  mole  of  A1  is  about  27  g,  so  0.3  kg  is  about  1 1  moles.  Using 
the  heat  of  fusion  from  Table  1,  it  takes  about  120  kJ  to  melt  this  size  sample  of 
aluminum.  Adding  this  number  to  the  previously  calculated  Q ,  it  takes  roughly  300  kJ  to 
melt  100  cm  of  aluminum. 

To  melt  a  sample  of  aluminum,  a  laser  must  impart  300  kJ,  and  do  so  at  a  rate 
greater  than  200  kW  in  order  to  exceed  the  cooling  by  radiation  and  conduction.  A  0.5 
MW  laser  would  accomplish  this  in  about  1  second.  Leaving  room  for  atmospheric 
losses  as  the  beam  propagates  to  the  target,  a  FEL  needs  to  be  on  the  order  of  1  MW 
power  level  to  be  an  effective  weapon. 
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III.  FREE  ELECTRON  LASER  COMPONENTS 


FELs  can  be  divided  into  two  main  categories  based  on  their  design;  oscillators 
and  amplifiers.  As  shown  in  Figure  2,  most  major  components  are  common  to  both 
configurations.  In  the  oscillator  configuration,  an  optical  resonator  cavity  is  formed  by 
highly  reflective  mirrors,  similar  to  traditional  molecular  lasers.  This  allows  the  optical 
wave  to  undergo  many  passes  through  the  undulator,  gaining  power  and  inducing 
stimulated  emission.  In  the  amplifier  configuration,  there  is  no  optical  resonator.  Instead 
a  seed  laser  is  used  to  induce  stimulated  emission.  Amplifiers  typically  have  longer 
undulators  then  oscillators  because  the  optical  wave  only  has  one  pass  in  which  to  gain 
power.  A  third  type  of  FEL  is  the  self-amplified  spontaneous  emission  (SASE)  FEL, 
which  is  essentially  an  amplifier  with  no  seed  laser. 


Superconducting  Electron  Accelerator 


Liquid  Helium  Refrigerator 


FEL  Oscillator 

|  I  ♦  BBBBfit  I  | 

Undulators 

#ftAAAA/V\AAA/\AMWVW^ 

FEL  Amplifier 


Figure  2. 


Basic  components  of  a  Free  Electron  Laser,  from  [1] 
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A. 


ELECTRON  BEAM  COMPONENTS 


Before  discussing  the  components  related  the  electron  beam,  it  is  instructive  to 
discuss  the  quality  of  the  electron  beam  itself.  An  important  measure  of  electron  beam 
quality  is  transverse  emittance,  which  is  the  product  of  the  RMS  width  and  the  angular 
spread  of  the  beam,  as  measured  at  the  beam  waist.  Electrons  traveling  outside  of  the 
optical  beam  will  not  contribute  to  the  optical  gain,  so  it  is  desirable  to  keep  the  electron 
beam  radius  smaller  than  the  waist  radius  of  the  optical  beam.  Another  important 
measure  of  electron  beam  quality  is  energy  spread.  It  is  desirable  for  all  of  the  electrons 
in  a  pulse  to  have  very  close  to  the  same  energy.  As  will  be  shown  in  chapter  IV,  the 
energy  of  the  electrons  entering  the  undulator  is  of  critical  importance  to  the  operation  of 
the  laser.  These  criteria  must  be  taken  into  account  when  considering  each  of  the 
following  components  in  the  electron  beam’s  path. 

1.  Injector 

The  injector  is  where  the  electron  beam  originates.  There  are  several  ways  of 
producing  free  electrons,  the  most  common  being  emission  from  a  photocathode  via  the 
photoelectric  effect.  A  drive  laser  is  used  to  eject  electrons  from  the  cathode  surface. 
This  drive  laser  is  usually  a  solid-state  pulsed  laser,  with  a  pulse  rate  around  700  MHz. 
The  ejected  electrons  enter  an  RF  injector  that  increases  their  energy  to  about  5  million 
electron  volts  (MeV).  This  type  of  injector  is  well  suited  for  high  power  FELs  because  it 
can  produce  short  electron  pulses  picoseconds  in  length,  and  can  produce  a  high  current 
electron  beam  of  excellent  quality. 

2.  Accelerator 

The  electron  beam  from  the  injector  goes  directly  into  the  accelerator  portion  of 
the  laser,  where  RF  energy  increases  the  electron  energy  from  around  5  MeV  to  around 
100  MeV.  The  accelerator  consists  of  several  cavities  made  of  pure  niobium,  which  at 
temperatures  around  2K,  becomes  superconducting.  The  alternating  RF  fields  of  the 
cavities  are  generated  by  a  klystron,  and  operate  at  the  same  frequency  as  the  electron 
beam  injector.  The  electron  pulses  coming  from  the  injector  are  in  phase  with  the  RF 
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fields,  so  they  accelerate  as  they  pass  through  the  cavities.  This  can  also  have  the  effect 
of  shortening  the  electron  pulses,  since  the  slower  electrons  can  be  accelerated  more  than 
the  faster  electrons. 

After  the  electrons  pass  through  the  undulator,  they  again  pass  through  the 
accelerator,  this  time  -180  degrees  out  of  phase  with  the  RF  fields.  This  causes  the 
electrons  to  give  up  most  of  their  energy  before  going  to  the  beam  dump.  This  second 
pass  through  the  accelerator  reduces  cooling  and  shielding  requirements  because  lower 
energy  electrons  produce  less  heat  and  radiation  upon  disposal.  Additionally,  the  energy 
given  up  by  the  electrons  during  their  second  pass  is  effectively  recycled,  lowering  the 
RF  input  power  required,  and  increasing  the  overall  efficiency  of  the  system. 

3.  Undulator 


The  undulator  is  contained  within  the  optical  resonator  cavity,  and  consists  of  a 
series  of  permanent  magnets  that  set  up  a  magnetic  field  of  alternating  polarity.  This 
alternating  field  is  perpendicular  to  the  direction  of  electron  motion,  and  induces  a  slight 
wiggle  to  the  electron’s  trajectory.  For  this  reason  the  undulator  is  sometimes  called  the 
“wiggler”.  As  the  relativistic  electrons  oscillate,  they  emit  synchrotron  radiation  along 
the  axis  of  the  undulator.  This  emitted  radiation  amplifies  the  optical  beam  already 
within  the  resonator  cavity.  Figure  3  depicts  this  process. 


Figure  3.  FEL  Undulator,  from  [23] 
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Undulators  may  have  either  linear  polarization,  as  shown  in  Figure  3,  or  helical 
polarization.  Linear  undulators  impart  a  sinusoidal  motion  to  the  electrons,  whereas 
helical  undulators  impart  a  corkscrew  motion.  In  either  case,  the  spacing  between  the 
magnets  is  called  the  undulator  period.  The  undulator  period  plays  a  role  in  determining 
the  frequency  of  the  electron  oscillations,  and  is  important  in  detennining  the  wavelength 
of  the  laser. 

4.  Beam  Dump 

Since  the  electron  beam  quality,  which  is  vital  to  the  operation  of  the  laser, 
degrades  as  the  electrons  pass  through  the  undulator  and  the  turning  magnets,  the  electron 
beam  cannot  be  used  in  the  undulator  a  second  time.  The  beam  dump  is  where  the 
electrons  end  their  journey,  after  passing  through  the  accelerator  a  second  time.  The 
second  pass  through  the  accelerator,  this  time  out  of  phase  with  the  RF  field,  reduces  the 
electron  energy  to  less  than  10  MeV. 

The  beam  dump  is  essentially  a  piece  of  metal,  usually  aluminum,  copper,  or 
tungsten  that  absorbs  the  electrons.  The  beam  dump  is  water  cooled  to  dissipate  the  heat 
resulting  from  the  absorption  of  the  electrons. 

5.  Electron  Beam  Transport  System 

Though  not  pictured  in  Figure  2,  two  types  of  magnets  may  be  found  throughout 
the  path  of  the  electron  beam.  Turning  magnets  are  used  to  redirect  the  electron  beam  so 
the  components  discussed  above  need  not  all  be  in  a  straight  line.  This  allows  the  whole 
system  to  be  more  compact.  They  are  sometimes  referred  to  as  bending  magnets  or 
dipole  magnets. 

Quadrupole  magnets  are  used  to  focus  the  electron  beam  and  maintain  high  beam 
quality,  as  discussed  at  the  beginning  of  this  section.  There  are  two  types  of  quadrupole 
magnets.  “F”  magnets  focus  the  beam  horizontally  while  defocus ing  it  vertically,  while 
“D”  magnets  have  the  opposite  effect.  How  these  magnets  work  will  be  discussed  in 
more  detail  in  Chapter  V. 
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The  path  of  the  electron  beam  is  kept  at  a  vacuum  to  prevent  the  electrons  from 
interacting  with  gas  molecules,  which  would  cause  an  unacceptable  reduction  in  beam 
energy  and  quality. 

B.  OPTICAL  BEAM  COMPONENTS 

1.  Optical  Resonator  Cavity 

As  mentioned  previously,  the  presence  of  an  optical  resonator  is  the  primary 
difference  between  oscillator  type  FELs  and  amplifier  type  FELs.  In  an  amplifier,  a  seed 
laser  synchronized  with  the  pulses  of  the  electron  beam  is  used  to  initiate  the  stimulated 
emission  process.  The  optical  beam  undergoes  significant  amplification  through  the 
entire  length  of  the  undulator,  but  makes  only  one  pass  through.  Therefore  an  amplifier 
FEL  has  no  optical  resonator. 

In  an  oscillator,  there  is  no  seed  laser.  The  optical  wave  is  initiated  by 
spontaneous  emission  from  the  electrons  into  the  resonator.  It  then  oscillates  back  and 
forth  through  the  undulator,  acting  as  the  seed  for  stimulated  emission.  The  mirror  on 
one  end  of  the  optical  cavity  is  partially  transparent  in  order  to  couple  the  laser  light  out 
of  the  cavity. 

Since  the  optical  beam  makes  many  passes  through  the  undulator,  the  oscillator 
configuration  can  use  a  shorter  undulator  than  the  amplifier  configuration.  This  means 
the  overall  system  can  be  more  compact.  This  is  not  without  trade-offs,  however.  First, 
the  distance  between  the  mirrors  of  the  optical  cavity  must  be  precisely  tuned  and 
maintained  so  the  successive  passes  of  light  pulses  remain  in  phase  with  the  electron 
pulses.  Second,  if  the  out-coupling  mirror  is  50%  transparent,  then  for  a  megawatt  class 
FEL,  the  optical  beam  within  the  resonator  cavity  will  be  several  megawatts.  If  this 
amount  of  power  is  confined  within  too  small  of  a  spot  size,  then  it  may  destroy  the 
mirrors  of  the  cavity  before  it  could  ever  be  directed  toward  a  target. 
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2. 


Beam  Director 


For  an  FEL  weapon  system,  the  beam  director  is  where  the  laser  light  actually 
leaves  the  ship  to  travel  to  the  target.  The  beam  director  must  be  able  to  aim  the  beam 
and  track  a  target.  This  capability  has  already  been  demonstrated,  along  with  the  use  of 
adaptive  optics  to  correct  for  atmospheric  turbulence  [24], 

3.  Optical  Beam  Transport  System 

The  optical  beam  must  be  transported  from  the  optical  cavity  to  the  beam  director. 
This  is  accomplished  with  a  series  of  mirrors,  and  under  vacuum  to  prevent  the  beam 
from  losing  energy  to  scattering  and  absorption.  The  beam  must  also  be  sufficiently  large 
so  as  not  to  destroy  the  mirrors  as  it  makes  its  way  to  the  beam  director,  similar  to  the 
problem  faced  within  the  optical  resonator. 

C.  AUXILIARY  SYSTEMS 

1.  Refrigeration  System 

As  mentioned  previously,  the  accelerator  must  be  kept  around  2  K  to  maintain 
superconductivity.  This  can  be  accomplished  with  a  combination  of  liquid  nitrogen  and 
liquid  helium.  This  cooling  system  can  also  provide  cooling  for  the  injector,  the  beam 
dump,  and  for  the  mirrors  of  the  optical  resonator  cavity. 

2.  Shielding 

The  injector,  undulator,  and  beam  dump  all  generate  radiation  from  which 
humans  need  to  be  shielded.  Because  the  electrons  are  recirculated  through  the 
accelerator  a  second  time,  their  energy  level  is  below  the  threshold  for  neutron  production 
(-10  MeV)  when  they  reach  the  beam  dump.  This  greatly  reduces  the  shielding 
requirements.  However,  the  electrons  entering  the  beam  dump  do  generate 
bremsstrahlung  radiation.  Additionally,  synchrotron  radiation  is  generated  at  the  bending 
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magnets  where  the  electron  beam  direction  is  altered.  The  radioactivity  decays  very 
rapidly  once  the  system  is  shut  down,  but  some  shielding  is  required  to  protect  personnel 
from  radiation  hazards  during  laser  operation. 

3.  Vibration  Control 

As  mentioned  previously,  mirror  spacing  and  alignment  is  crucial  to  FEL 
operation.  Vibrations  normally  found  in  a  shipboard  environment  could  interfere  not 
only  with  mirror  alignment,  but  also  with  the  alignment  of  the  quadrupole  magnets. 
Vibration  isolation  mounting,  and  an  active  alignment  system  will  be  required  to  ensure 
proper  operation. 
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IV.  BASIC  PHYSICS  OF  FREE  ELECTRON  LASER  OPERATION 


A.  ELECTRON  BEHAVIOR 

In  the  following  sections,  the  z-axis  will  denote  the  longitudinal  axis,  passing 
through  the  center  of  the  undulator.  Electron  velocities  will  usually  be  normalized  to  the 
speed  of  light  as  follows:  Equation  Chapter  4  Section  1 

Px  =  VJC 

Py  =  Vy/C 

P,  =  vJc 

P=(P^P^P-) 

The  undulator  period  is  20 . 

1.  Undulator  Fields  and  Electron  Motion 

To  first  order  approximation,  electrons  travel  straight  through  the  undulator,  down 
the  z-axis  at  very  near  the  speed  of  light 

z(t )  =  v.t  =  fi_ct  «  ct.  (4. 1) 

However,  there  is  a  more  subtle  motion  undergone  by  the  electrons  due  to  the  fields  that 
are  present.  Inside  the  undulator,  the  electrons  are  exposed  to  three  fields;  the  magnetic 
field  from  the  undulator  magnets,  and  the  magnetic  and  electric  fields  of  the  optical  wave. 
The  magnetic  field  from  the  magnets  of  a  helical  undulator  is 

Bm  =  B(cos(k0z),  sin(A:0z),0),  (4.2) 

where  k0  =2  is  the  undulator  wave  number.  The  electric  and  magnetic  fields  from 
a  corresponding  circularly  polarized  optical  plane  wave  are 

E0  =  E  (cos  i//,- sin  i//, 0^,  (4.3) 

B0=  E  (sin  i//,  cos  if/,  0),  (4.4) 

where  y/  =  kz  -  cot  +  (f> ,  co  is  the  optical  frequency,  (f)  is  the  optical  phase,  and  k  =  2k/ X 
is  the  optical  wave  number. 
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Electron  motion  in  the  presence  of  these  fields  can  be  detennined  using  the 
relativistic  Lorentz  force  equations, 


djyP) 

dt 


dy 

dt 


7  = 


—  (E  +  fixB), 
mec 


e 

rnec 


P-E, 


(4.5) 

(4.6) 

(4.7) 


Where  e  is  the  electron  charge  magnitude,  me  is  the  electron  mass,  and  y  is  the  Lorentz 

factor.  Substituting  the  field  equations  (4.2,  4.3,  and  4.4)  into  the  Lorentz  force  equation 
(4.5)  results  in  the  following  equations  of  electron  motion, 

=  -^[Ecos  i/(l  -  fiz)  -  PZB  sin(V)] ,  (4.8) 

dt  mec 
d(yj3)  e 

— —  = - [~E  sin  i/(l  -  J3Z)  +  pzB  cos  (k0z)\ ,  (4.9) 

dt  mec 

=  — —\E(px  cos (i//)  -  J3V  sin(i/))  +  B(J3X  sin (k0z)  -  J3V  cos (k0z)\ .  (4.10) 

dt  mec 

Since  the  electrons  are  traveling  close  to  the  speed  of  light,  /?  «  1,  so  that  ( 1  -  fi. ) 


is  very  small.  Therefore  the  first  term  inside  the  brackets  in  (4.8)  and  (4.9)  can  be 
neglected.  Physically,  this  means  the  electrons  are  traveling  so  fast  that  the  influence  of 
the  optical  wave’s  electric  field  on  electron  transverse  trajectories  is  negligible  compared 
to  the  influence  of  the  undulator.  Equations  (4.8)  and  (4.9)  can  then  be  integrated  with 
respect  to  time  to  produce 

Px  =~— cos(k0z),  (4.11) 

Y 

Py  =  -  —  sin(k0z),  (4-i2) 

Y 

where  K  is  the  dimensionless  undulator  parameter,  defined  as  K  =  eB T0/ 2  7imec2  .  The 
constants  of  integration  have  been  taken  to  be  zero,  corresponding  to  perfect  injection  of 
the  electrons  into  the  helical  undulator.  Noting  that  k0z  «  k0ct  =  co0t,  (4.1 1)  and  (4.12) 
can  be  integrated  to  obtain 
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(4.13) 


x(t )  = - sin(<u0t)  = - -sin(<u0t), 

a>0y  2  ny 

Kc  KA, 

y(t)  =  — -cos(ny)  =  — 5 -cos(ny),  (4.14) 

a>0y  2  ny 

which  gives  the  trajectory  of  the  electrons  in  the  x  and  y  directions.  In  a  helical 
undulator,  the  motion  of  the  electrons  resembles  a  corkscrew.  It  is  instructive  to  use 
some  typical  values  to  put  this  corkscrew  motion  into  perspective.  Typically,  K  ~  1, 
y  ~  100,  and  A0  ~  5  cm.  This  gives  oscillations  in  the  x  and  y  directions  on  the  order  of 
80  microns,  which  is  almost  half  of  a  typical  electron  beam  waist  diameter  of  200 
microns.  The  frequency  of  these  oscillations  is  <x>0/2tt  =  c/A0  =6  GHz.  The  electrons 
are  oscillating  at  microwave  frequencies,  but  the  laser  operates  at  infrared  frequencies 
due  to  relativistic  Doppler  shifts  [1]. 


2.  Microscopic  Electron  Motion  and  the  Pendulum  Equation 


The  longitudinal  motion  of  the  electrons  is  very  straightforward,  as  described  by 
(4.1).  In  the  previous  section,  the  transverse  motion  of  the  electrons  was  examined  on  a 
scale  where  the  effects  of  the  optical  wave  fields  were  considered  negligible.  Now  we 
will  investigate  electron  motion  on  an  even  smaller  scale,  where  it  will  be  shown  that  the 
effects  of  the  optical  wave  fields,  while  small  in  magnitude,  are  vital  to  understanding 
FEL  operation. 


Substituting  (4. 1 1)  and  (4. 12)  into  (4. 10)  gives  us 

K  ,  K 


d(ypz) 

dt 


m  c 


+B 


- cos  k0z  cos  y/  -\ - sin  k0x  sin  yr 

v  r  r  j 

-  —  cos  k0z  sin  k0z  +  —  sin  k0z  cos  k0z 


Y 


r 


j 


eKE 


y/3z  +y/T  = - (cosk0zcos^r  +  sinA:0xsin^'), 

ymec 

YP,  =^Lcos(k0z  +  y/)~ yfiz . 
ymec 


We  can  also  substitute  (4.1 1)  and  (4.12)  into  (4.6)  to  get 


(4.15) 


(4.16) 

(4.17) 
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(4.18) 


dy  .  eKE  \ 

—  =  y  = - cos  ( k0z  +  y/  )  ■ 

dt  ymec 


If  (4.18)  is  then  substituted  into  (4.17),  we  get 


A=(i-A) 


eKE 

y2mec 


cos 


[k0z  +  y/). 


(4.19) 


At  this  point  it  is  useful  to  introduce  the  electron  phase,  £  =  (k  +  k0)z  -  cot .  Since 
k ,  k0,  and  co  are  fixed,  £  follows  the  electron’s  position  in  z  .  The  electron  phase  may 
be  interpreted  as  a  measure  of  the  electron’s  position  relative  to  the  optical  wave.  It  is 
this  relative  position  that  determines  whether  the  electron  gives  energy  to  the  optical 
wave,  or  vice  versa.  Taking  the  first  and  second  derivatives  with  respect  to  time  of  the 
electron  phase,  we  get 

£=(k  +  k,)P:c-a>,  (4.20) 

C  =  (k  +  ko)c0z.  (4.21) 


Solving  for  /?,  and  jdz  we  get 


A  = 
A  = 


£  +  C0 
(. k  +  k0)c 

c 


(4.22) 

(4.23) 


(. k  +  k0)c 

Since  A0,  it  follows  that  k  k0 .  Also,  since  fd_  «  1  and  co  =  kc ,  we  see  from  (4.22) 
that  £  co .  So  (4.22)  can  then  be  simplified  to 

co 


A  = 


Substituting  (4.23)  and  (4.24)  into  (4.19) 


(k  +  k0)c 


(4.24) 


C 


(k  +  k0)c 


1- 


co 


(k  +  k0)c  )y2mec 


eKE 


cos 


{C  +  <?>)’ 


f=  eKE 
y  me 


\_{k  +  K)  c-coj  cos(^  +  (f>)  , 


eKEk„ 

y2me 


-cos 


{£+*)■ 


(4.25) 
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Before  proceeding,  it  will  be  useful  to  nonnalize  our  parameters  to  the  dimensions 
of  the  undulator.  If  there  are  N  periods  in  the  undulator,  then  the  undulator  length  is 
L  =  N/ 10 ,  and  we  can  introduce  the  dimensionless  time  parameter  r  =  ct  /  L  .  As  the 
optical  wave  makes  a  single  pass  through  the  undulator,  z  goes  from  0  to  1 .  Derivatives 

O 

may  now  be  taken  with  respect  to  r  rather  than  t ,  and  denoted  as  (..)  =  d(..)/dz  rather 
than  (..)  =  d(..)/dt . 


It  is  also  convenient  to  combine  the  constants  in  front  of  the  cosine  term  in  (4.25) 
to  introduce  the  dimensionless  optical  field  magnitude 

,  ,  InNeKEL 


\a\  = 


2  2 

Y  mec 


With  these  simplifications,  (4.25)  may  now  be  written  as 

00  o 

£=v  =  |a|cos(<^+  (f>), 

O 

where  v  =  C,  is  the  electron  phase  velocity. 


(4.26) 


(4.27) 


Equation  (4.27)  is  in  the  form  of  the  well-known  pendulum  equation,  and 
describes  the  microscopic  behavior  of  the  electrons  within  an  optical  wavelength.  The 
behavior  of  simple  pendulums  is  well  understood,  and  can  be  easily  explained  with  the 
use  of  phase  space  diagrams.  Likewise,  (4.27)  allows  the  evolution  of  electrons  within 
the  undulator  to  be  demonstrated  with  phase  space  diagrams. 


3.  Phase  Space  Evolution,  Resonance,  and  Bunching 

The  meaning  of  electron  phase  was  described  above.  Equation  (4.27)  describes 
how  the  electron  phase  velocity  changes  with  time.  Changes  in  electron  phase  velocity 
result  from  exchanges  of  energy  between  the  electrons  and  the  optical  wave  [25].  The 
electron  phase  velocity  increases  when  the  electron  gains  energy,  and  decreases  when  the 
electron  loses  energy.  If  there  is  no  energy  exchange,  the  electron  phase  remains  constant 
and  the  electron  phase  velocity  is  zero.  The  mechanical  analog  of  this  is  a  motionless 
pendulum  at  the  bottom  of  its  arc,  with  no  exchange  between  kinetic  and  potential 
energy.  Recalling  that  <x>  =  kc ,  we  can  set  (4.20)  equal  to  zero  and  solve  to  find 
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which  is  known  as  the  resonance  condition.  Note  that  the  optical  wavelength,  A,  depends 
upon  the  undulator  wavelength,  A0.  This  is  one  of  the  features  that  make  an  FEL 
“tunable.”  If  we  replace  //.  with  v./c  in  (4.28),  then  we  see 


;  l~vJc  -  i  c~v-- 

sL  An  ,  /Lr\ 


vjc 


=  T0(c/v,-l). 


A  +  Ag  Ag 


(4.29) 


The  right  hand  side  of  (4.29)  is  simply  the  amount  of  time  required  for  an  electron  to  pass 
through  one  undulator  wavelength.  At  resonance,  this  is  equal  to  the  amount  of  time 
required  for  the  optical  wave  to  travel  one  optical  wavelength  further.  In  other  words,  at 
resonance  the  optical  wave  outruns  the  electron  by  a  distance  of  one  optical  wavelength. 
Figure  4  illustrates  this  concept.  The  red  circle  represents  a  single  electron  as  it  travels 
from  left  to  right.  The  undulator  wavelength  is  shown  in  green,  and  the  optical  wave  is 
shown  in  blue. 


Figure  4.  Electron  in  resonance  condition,  from  [1] 

Figure  5  depicts  the  evolution  in  phase  space  of  20  sample  electrons  as  they  pass 
through  the  undulator  from  left  to  right.  The  electrons,  shown  as  red  dots,  begin  with 
zero  initial  phase  velocity,  in  other  words  on  resonance,  but  evenly  distributed  in  phase. 
In  the  phase  space  diagrams,  the  vertical  axis  is  phase  velocity,  which  corresponds  to 


electron  energy.  The  horizontal  axis  is  electron  phase,  which  represents  electron  position 
relative  to  the  optical  wave.  The  “lemon-shaped”  curve  represents  the  separatrix,  which 
carries  the  same  meaning  here  as  it  would  in  a  phase  space  plot  of  pendulum  motion,  with 
paths  inside  representing  closed  orbits  and  paths  outside  representing  open  orbits.  In  this 
case,  all  electrons  are  following  closed  orbits.  Note  that  at  the  end  of  the  undulator,  half 
of  the  electrons  have  decreased  in  phase  velocity,  meaning  energy  was  transferred  to  the 
optical  wave,  and  half  have  increased  in  phase  velocity,  meaning  energy  was  taken  away 
from  the  optical  wave.  At  resonance  there  is  no  net  optical  gain. 


•lmctron  beam 
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Figure  5.  Phase  space  evolution  at  resonance,  from  [1]. 

Also  worth  noting  is  the  bunching  effect.  When  entering  the  undulator,  the 
electrons  are  uniformly  distributed  in  phase.  By  the  end  of  the  undulator,  the  electrons 
have  bunched  around  a  specific  value  of  phase.  The  reason  for  this  can  be  understood  by 
examining  (4.27).  Electrons  with  phases  in  the  range  -nj2  <  £+  <j)  <  n/2  will  yield  a 
positive  cosine  tenn,  resulting  in  a  positive  phase  acceleration,  and  an  increase  in  phase 
velocity.  Electrons  with  phases  in  the  range  n/2  <£+(j)<  In  12  will  yield  a  negative 
cosine  term,  resulting  in  a  negative  phase  acceleration,  and  a  decrease  in  phase  velocity. 
At  the  resonance  condition,  electrons  tend  to  bunch  around  the  value  of  phase 
corresponding  to  zero  phase  velocity. 
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A  more  favorable  result  in  terms  of  laser  operation  is  found  by  operating  above 
resonance,  which  can  be  achieved  by  injecting  electrons  with  slightly  higher  energy 
(~1%)  then  would  be  dictated  by  the  resonance  condition.  This  case  is  depicted  by 
Figure  6.  Note  that  in  this  example,  by  the  end  of  the  undulator  there  has  been  more 
downward  motion  in  phase  space  then  upward  motion.  This  means  that  overall,  the 
electron  pulse  has  transferred  energy  to  the  optical  wave,  resulting  in  a  net  optical  gain. 
Also  note  that 


•l«Ctron  beam 


c 


Figure  6.  Phase  space  evolution  above  resonance,  from  [1]. 

the  electrons  are  bunching  around  a  larger  phase  value  then  they  did  at  the  resonance 
condition.  In  terms  of  optical  gain,  it  is  better  to  bunch  electrons  around  a  phase  value 
that  causes  them  to  give  up  energy. 

Figure  7  illustrates  the  idea  of  bunching  electrons  in  a  region  that  is  favorable  to 
optical  gain.  The  red  region  represents  values  of  phase  corresponding  to  energy  transfer 
from  the  optical  wave  to  the  electrons,  or  absorption.  The  green  region  represents  values 
of  phase  corresponding  to  energy  transfer  from  the  electrons  to  the  optical  wave,  or  gain. 
The  upper  picture  shows  electrons  at  resonance,  as  in  Figure  5.  They  begin  at  ti 
uniformly  distributed  in  phase,  and  end  at  t2  bunched  around  the  phase  corresponding  to 
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zero  net  gain.  The  lower  picture  shows  electrons  above  resonance,  as  in  Figure  6.  They 
begin  at  ti  uniformly  distributed  in  phase,  and  end  at  t2  bunched  around  a  higher  value  of 
phase,  corresponding  to  positive  net  gain. 


Figure  7. 


Proper  electron  bunching  results  in  optical  gain,  from  [23]. 


B.  OPTICAL  WAVE  EQUATION 

The  dimensionless  optical  field  magnitude  was  introduced  in  (4.26).  The 
corresponding  wave  equation  can  be  expressed  as 

a  =  -  j  <  e~l(  > ,  (4.30) 

where  j  =  8N(e7rKL)2 p/y3  mec2  is  the  dimensionless  electron  beam  current,  p  is  the 
electron  beam  particle  density,  and  < ..  >  denotes  an  average  over  all  the  electrons  [1]. 
Note  that  if  the  electrons  are  evenly  distributed  throughout  phase,  the  average  is  zero. 
This  average  will  only  be  nonzero  when  the  electrons  are  bunched  around  a  phase  that 
has  a  non-zero  phase  velocity.  Therefore,  the  amount  by  which  the  optical  field  changes 
depends  on  the  electron  current,  and  the  electron  bunching. 
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1. 


GAIN  AND  SATURATION 


The  optical  beam  gain  is  the  fractional  increase  of  optical  field  energy  during  one 
pass  through  the  undulator,  defined  by 


a(  r) 

,|2 \t 

r(0)2| 

a(0): 

>  | 

(4.31) 


From  (4.27),  it  is  clear  that  a  strong  optical  field  (\a\>  7t)  is  necessary  for  good  electron 


bunching.  Provided  we  have  bunching  at  a  favorable  phase,  (4.30)  implies  that  the 
optical  field  will  grow,  leading  to  more  bunching.  The  interdependence  of  (4.27)  and 
(4.30)  constitutes  a  feedback  cycle  that  results  in  exponential  gain. 


If  the  optical  field  is  too  strong,  however,  the  average  electron  phase  can  evolve  to 
a  point  where  the  phase  velocity  begins  to  increase,  which  takes  energy  away  from  the 
optical  field.  This  causes  the  gain  to  decrease,  and  the  optical  growth  process  becomes 
saturated.  Figure  8  illustrates  what  this  “overbunching”  looks  like  in  phase  space. 


V 

S 

Figure  8.  Overbunching,  from  [1] 
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2.  Rayleigh  Length  Considerations 


In  order  for  a  FEL  to  have  gain,  there  must  be  overlap  between  the  electron  beam 
and  the  optical  beam.  To  quantify  this  overlap,  let  us  define  the  filling  factor  as  the  ratio 
of  electron  beam  area  to  optical  beam  area. 


Rayleigh  length,  Z0 ,  sometimes  called  depth  of  focus,  is  the  distance  from  the 

beam  waist  where  the  beam  radius  has  increased  by  a  factor  of  V2  [26],  For  a  circular 
beam,  this  equates  to  doubling  the  cross  sectional  area  of  the  beam.  Clearly,  the  filling 
factor  will  depend  upon  the  Rayleigh  length. 


If  we  assume  the  optical  beam  consists  only  of  the  fundamental  Gaussian  mode, 
then  it  can  be  shown  by  averaging  the  filling  factor  over  the  length  of  the  undulator  that 
the  gain  is  [27] 


Go  A/3 

Z0/L  +  L/(\2Z0y 


(4.32) 


where  G0  is  the  maximum  gain  at  the  optimal  Rayleigh  length.  This  expression  shows 
that  the  gain  will  go  down  as  the  Rayleigh  length  gets  shorter. 


The  shorter  the  Rayleigh  length,  the  more  divergent  is  the  optical  beam.  If  the 
optical  beam  is  not  sufficiently  divergent,  then  there  may  be  too  much  power 
concentrated  in  too  small  an  area  to  be  handled  by  the  mirrors  of  the  resonator  cavity.  As 
will  be  shown  later,  Rayleigh  length  also  affects  the  laser’s  sensitivity  to  changes  in  the 
electron  beam’s  trajectory. 
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V.  ELECTRON  BEAM  OPTICS 


In  chapter  III,  the  electron  beam  transport  system  was  mentioned  along  with  the 
types  of  magnets  commonly  used  in  FELs.  In  this  chapter  we  will  examine  in  detail  how 
these  magnets  affect  the  trajectory  of  an  electron.  In  our  calculations  we  will  use  a 
coordinate  system  that  follows  the  electron.  The  direction  of  electron  travel  will  be 
designated  s  ,  while  the  horizontal  and  vertical  planes  are  defined  as  x-s  and  y-s 
respectively.  The  right-hand  side  of  Figure  9  illustrates  this  coordinate  system,  where  s 
comes  out  of  the  page.  Equation  Chapter  5  Section  1 

A.  QUADRUPOLES 

As  previously  stated,  quadrupoles  are  used  to  focus  an  electron  beam.  Figure  9 
shows  a  cross  section  of  an  “F”  type  quadrupole,  which  is  used  to  focus  an  electron  beam 
in  the  horizontal  plane.  “D”  type  quadrupoles,  used  to  focus  an  electron  beam  in  the 
vertical  plane,  are  identical  to  the  “F”  type  except  that  the  poles  are  reversed. 


> 


Figure  9.  Cross-section  of  an  F  type  quadrupole  magnet,  from  [28]. 

In  calculating  the  effect  quadrupoles  have  on  an  electron’s  trajectory,  it  will  be  assumed 
that  the  magnetic  field  is  constant  along  the  s-axis,  and  zero  outside  the  magnet.  Results 
obtained  through  these  assumptions  have  been  shown  to  be  in  excellent  agreement  with 
measurements  [5]. 
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The  Lorentz  forces  experienced  by  an  electron  located  at  (x,  v) ,  traveling  in  the  s 


d2x 

B  e 

(5.1) 

dt 2  ~ 

xvs 

a 

d2y 

B\e\ 

(5.2) 

y  _ 

e  dt 2  " 

yvs 

a 

direction  with  velocity  v, ,  and  passing  through  an  “F”  type  quadrupole  are 

Fx  =  yme 

Fy  =  Yme 

where  B  is  the  magnetic  field  magnitude  at  the  pole  face  and  a  is  the  semi-aperture,  as 
shown  in  figure  9.  Note  that  if  x  is  positive,  the  electron  experiences  a  force  in  the  -x 
direction,  and  vice  versa.  This  means  the  quadrupole  focuses  in  the  horizontal  plane. 
Also  note  that  the  quadrupole  is  defocusing  in  the  vertical  plane.  If  the  electron  is  in  the 
center  of  the  quadrupole,  no  force  is  experienced. 


It  will  be  convenient  to  replace  the  time  derivatives  in  (5.1)  and  (5.2)  with 
derivatives  with  respect  to  s .  This  is  possible  because  the  electron’s  position  in  s  is 
uniquely  defined  for  any  given  time  by 

s  =  vst  ~  ct , 

therefore 


d 1  2  d 1 


dt 


,  =  c  , 
ds 


Substituting  (5.3)  into  (5.1)  and  (5.2)  results  in 

d2x 


ds 
d2y 
ds2 


=  ~KX  , 


Ky 


(5.3) 

(5.4) 

(5.5) 


where  k  =  \e\ B/ymeac  is  the  normalized  linear  field  gradient.  Equations  (5.4)  and  (5.5) 


are  second-order  linear  homogeneous  differential  equations  that  are  easily  solved 
analytically.  Taking  the  initial  conditions  at  5  =  0  to  be  x  =  x0 ,  y  =  y0 ,  dx/ds  =  x'  =  x'0, 
and  y'  =  y'0 ,  the  solutions  are 


x  =  x0  cos 


+ 


y  = 


(5.6) 

(5.7) 
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x  =— X, 


ly[icsm^Sy[fc'j  +  x’0  cos^sJk  j  , 


y'  =  y0  V/c  sinh  (s4k  y'0  cosh  (vV/c  ). 

A  more  convenient  way  of  expressing  these  equations  is  in  matrix  notation 


(5.8) 

(5.9) 


Vx7 


(y\ 

vV'y 


r 


cos 


(syfic^ 


-^=  sin  (s  V/c 'j 

yJtC  V  ' 

IVv 

cos 

IxJ 

) 

^=sinh(sV/c 

yj/C  v 

)  rv 

cosh^V/r) 

w 

) 

(5.10) 


(5.11) 


If  s  in  (5.10)  and  (5.11)  is  taken  to  be  the  effective  length  of  the  quadrupole,  then  these 
matrices  enable  us  to  calculate  the  final  position  and  trajectory  of  an  electron  whose 
initial  position  and  trajectory  are  known  [5],  For  this  reason,  they  are  appropriately 
tenned  transfer  matrices,  analogous  to  what  may  be  found  in  the  field  of  optics  to 
describe  lenses.  The  transfer  matrices  for  a  “D”  type  quadrupole  are  the  same,  simply 
switching  x  with  y  and  x'  with  /  . 


If  we  wish  to  investigate  the  effect  of  a  displaced  quadrupole,  we  need  only 
transfonn  the  electron  to  a  frame  where  the  magnet  is  centered,  apply  the  transfer 
matrices  defined  above,  then  transform  the  electron  back  to  the  original  reference  frame. 
As  an  example, 

x  =T-(x0  +  d)-d  ,  (5.12) 

where  T  is  a  transfer  matrix  and  d  is  the  quadrupole  displacement.  Displaced 
quadrupoles  alter  the  trajectory  of  the  electron  beam,  which  can  effect  the  operation  of 
the  FEL.  This  will  be  explored  in  more  detail  later. 


B.  DIPOLES 


Dipole  magnets  are  used  to  steer  the  electron  beam.  They  allow  the  beam  to  be 
doubled  back  so  the  overall  system  takes  up  less  space.  Since  their  only  purpose  is  to 
change  the  direction  of  the  electron  beam,  there  are  only  two  magnetic  poles,  yielding  a 
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magnetic  field  with  only  one  directional  component.  We  will  assume  this  is  in  the 
vertical  plane.  The  relativistic  Lorentz  force  experienced  by  an  electron  inside  a  dipole  is 

Fx  =  ymec2x"  =  -evsBy  .  (5.13) 

Note  that  this  describes  the  force  in  the  laboratory  frame  of  reference.  If  the  electron  is 
precisely  on  the  s-axis  with  no  displacement  in  x  or  y,  it  will  remain  precisely  on  the  s- 
axis  with  no  transverse  displacement  because  of  the  way  the  s-axis  is  defined.  We  can 
think  of  dipoles  as  changing  the  path  of  s,  the  nominal  path  of  the  electron,  within  the 
laboratory  reference  frame.  Equation  (5.13)  causes  the  s-axis  to  curve  in  the  horizontal 
plane.  The  radius  of  curvature  is  found  by  equating  the  Lorentz  force  to  a  centripetal 
force, 

FR  =  Yme  ^  =  ~eVsBy  ■  (5-14) 

Solving  for  the  radius  of  curvature,  R  ,  gives 


r>neVs 

eBy 


(5.15) 


If  the  electron  is  not  precisely  on  the  s-axis,  but  rather  has  a  displacement  in  x ,  the  there 
will  be  an  apparent  acceleration  in  x  relative  to  the  s-axis.  In  the  laboratory  frame,  the 
displaced  electron  will  have  the  same  radius  of  curvature  as  an  electron  with  no 
displacement.  However,  in  our  co-moving  coordinate  system,  the  displaced  electron  will 
move  relative  to  the  non-displaced  electron.  This  motion  can  be  characterized  by  framing 
the  problem  as  a  mismatch  between  the  Lorentz  force  and  a  centripetal  force  with 
r  =  R  +  x . 


ymec2x"  =  Fr-  Florentz , 

9  (  /  \ 

ymec  x  =yyme—  J  - [^evsBy ), 

(  v ■  '  (•  ■  «). 


ymec2x"  =  | ym 


(R  +  x) 

Inserting  (5.15)  into  (5.16)  and  solving  for  x"  yields 


+ 


x  =—- 


R2 

The  solution  to  (5.17)  expressed  in  matrix  form  is 


(5.16) 

(5.17) 
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(5.18) 


To  obtain  the  transfer  matrix  for  the  vertical  direction,  recall  that  Bx  =  0 ,  which 
means  \/R  =  0 .  Substituting  this  into  (5.18)  gives 


1 


(5.19) 


\y'J  'vO  IAjv 

A  transfer  matrix  of  this  form  describes  a  drift,  or  region  where  there  are  no  forces  acting 
on  the  electron  in  that  direction. 


C.  EDGES 

In  deriving  the  transfer  matrices  for  a  dipole  magnet,  it  was  once  again  assumed 
that  the  magnetic  field  was  constant  inside  the  magnet,  and  zero  outside  the  magnet,  just 
as  had  been  done  for  the  quadrupole.  Furthermore,  it  was  assumed  that  the  magnetic 
field  began  and  ended  at  the  same  value  of  s  over  the  entire  x-s  plane.  This  is  only  true 
if  the  faces  of  the  dipole  are  perpendicular  to  the  s  axis.  This  is  true  for  sector  magnets 
as  shown  in  Figure  10.  However,  rectangular  magnets  are  often  used  because  they  are 
easier  to  manufacture  [5].  As  Figure  10  shows,  the  face  of  a  rectangular  magnet  is  tilted 
through  the  angle  y/ ,  which  affects  the  electron’s  trajectory  as  it  enters  and  exits  the 
dipole. 


sector  magnet 


rectangle  magnet 


x 


Figure  10.  Sector  versus  rectangular  magnets,  from  [5] 
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In  the  right  hand  side  of  Figure  10,  we  see  that  an  electron  with  initial  horizontal 
displacement  x(}  travels  a  shorter  distance  inside  the  magnet  than  an  electron  with  zero 
initial  displacement.  The  difference  in  path  length  is 

A/ =  x0  tan  ^ ,  (5.20) 


which  causes  the  electron  to  be  bent  through  a  smaller  angle.  This  difference  is 

A  /  tan  y/ 

A  a  =  —  =  xn - . 

R  R 


(5.21) 


After  passing  through  an  edge,  the  final  horizontal  displacement  of  the  electron  is  the 
same  as  the  initial  displacement,  but  the  trajectory  has  been  altered  by  A  a  .  Hence,  the 
transfer  matrix  for  the  edge  of  a  rectangular  dipole  magnet  is 


(5.22) 


(x^ 

1 

f  1  O' 

(x  \ 

= 

tan  y/ 

Wj 

, - —  1 

\x'0) 

V  R  J 

This  implies  the  vertical  component  of  the  magnetic  field  has  a  gradient  in  the  x 
direction.  By  symmetry  (Vx2?  =  0),  the  horizontal  component  of  the  magnetic  field 
must  have  a  gradient  in  the  y  direction  that  is  equal  in  magnitude  but  opposite  in 
direction.  Therefore, 


(5.23) 
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VI.  SIMULATIONS  TO  EXPLORE  RAYLEIGH  LENGTH 
DEPENDENCE  OF  GAIN  VERSUS  TILT 


As  mentioned  in  chapter  IV,  material  considerations  motivate  us  to  develop  short 
Rayleigh  length  FELs.  Also  as  mentioned  in  chapter  IV,  decreasing  the  Rayleigh  length 
can  affect  the  on  gain.  Experiments  at  Thomas  Jefferson  National  Accelerator  Facility 
(TJNAF)  were  conducted  to  explore  FEL  performance  as  it  relates  to  shift  and  tilt 
tolerances,  and  the  affect  that  Rayleigh  length  has  on  these  tolerances.  In  this  section  we 
discuss  simulations  that  were  performed  for  the  same  purpose. 

To  clarify  our  tenninology,  shift  refers  to  a  translation  of  the  electron  beam  along 
either  the  vertical  or  the  horizontal  axis,  Ax  or  Ay  .  Tilt  refers  to  rotation  of  the  beam 
about  a  point,  Ax'  or  Ay' .  In  aviation  parlance,  a  vertical  tilt  corresponds  to  pitch,  while 
a  horizontal  tilt  corresponds  to  yaw.  In  this  section  all  shifts  and  tilts  are  in  the  vertical 
direction,  which  is  in  the  plane  of  the  undulator  field. 

A.  3D  FEL  OSCILLATOR  SIMULATIONS 

The  3D  FEL  simulator  was  created  by  Professor  Colson  and  Professor  Blau  at  the 
Naval  Postgraduate  School  (NPS)  [29].  This  code  runs  on  a  64  node,  128  processor 
Apple  Xserve  cluster.  An  input  file  specifies  parameters  for  the  electron  beam,  the 
undulator,  and  the  optical  cavity,  and  the  output  provides  information  for  the  evaluation 
of  laser  performance,  such  as  gain  and  extraction.  Gain  was  defined  in  chapter  IV. 
Extraction  is  the  output  optical  power  divided  by  the  average  power  of  the  electron  beam. 
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Table  2  summarizes  the  properties  of  the  electron  beam,  which  match  the 
properties  of  the  system  at  TJNAF.  The  energy  of  the  electron  beam  exiting  the 
accelerator  is  1 1 1  MeV.  The  electron  pulse  rate  is  4.7  MHz,  with  a  pulse  length  of  0.12 
mm. 


Eb 

ELECTRON  BEAM  PARAMETERS 

Beam  energy  (MeV) 

111 

qbunch 

Bunch  charge  (nC) 

0.108 

rbx 

Beam  radius,  x  rms  (mm) 

0.2 

rby 

Beam  radius,  y  rms  (mm) 

0.2 

tb 

Pulse  duration,  FWHM  (ps) 

0.4 

prf 

Pulse  rep  frequency  (MHz) 

4.7 

lb 

Pulse  length,  FWHM  (mm) 

.12 

gamma 

Lorentz  factor 

218 

Ipeak 

Peak  current  (A) 

270 

Iavg 

Average  current  (mA) 

0.5 

emitx 

Normalized  rms  x  emittance  (mm  mrad) 

7.5 

emity 

Normalized  rms  y  emittance  (mm  mrad) 

7.5 

emitl 

Longitudinal  emittance  (keV  ps) 

70 

dgog 

Beam  energy  spread  (%) 

0.37 

dthetax 

Beam  angular  spread,  x  rms  (mrad) 

0.17 

dthetay 

Beam  angular  spread,  y  rms  (mrad) 

0.17 

rho 

Beam  density  (l/cmA3) 

2.2xl013 

Pb 

Beam  average  power  (MW) 

0.46 

Table  2.  Electron  Beam  Parameters 


The  properties  of  the  undulator  are  summarized  in  Table  3,  which  again  were 
chosen  to  correspond  with  the  system  at  TJNAF.  There  are  30  undulator  periods,  with  an 
undulator  period  of  5.5  cm,  to  give  an  undulator  length  of  165  cm. 


lambdaO 

UNDULATOR  PARAMETERS 

Undulator  period  (cm) 

5.5 

N 

Number  of  periods 

30 

gap 

Undulator  gap  (cm) 

4 

Bpeak 

Undulator  peak  magnetic  field  (T) 

0.22 

Brms 

Undulator  magnetic  field,  rms  (T) 

0.16 

K 

Undulator  parameter,  rms 

0.81 

L 

Undulator  length  (cm) 

165 

Table  3.  Undulator  Parameters 
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Table  4  lists  the  optical  cavity  parameters.  The  optical  cavity  length,  which  is  the 
distance  between  the  mirrors,  is  32  m.  The  cavity  quality  factor  is  17,  giving  an  output 
coupling  of  5.9%.  The  optical  wavelength  is  0.96  microns. 


OPTICAL  CAVITY  PARAMETERS 


s 

Cavity  length  (m) 

32 

zo 

Rayleigh  length  (cm) 

70 

loss 

Mirror  losses  per  pass  (%) 

6 

lambda 

Optical  wavelength  (microns) 

0.96 

WO 

Mode  waist  radius,  1/e  (mm) 

0.5 

Wmir 

Mode  radius  at  mirrors,  1/e  (cm) 

1.1 

Qn 

Quality  factor 

17 

Table  4.  Optical  Cavity  Parameters 


Table  5  lists  the  dimensionless  parameters  that  are  the  direct  inputs  into  the 
simulation.  Their  values  are  based  on  the  physical  parameters  listed  in  Tables  2,  3,  and  4. 


DIMENSIONLESS  PARAMETERS 


i 

Normalized  current  density,  linear  undulator 

20 

sx 

Normalized  beam  radius 

0.28 

sy 

Normalized  beam  radius 

0.28 

stx 

Normalized  beam  angular  spread 

0.40 

sty 

Normalized  beam  angular  spread 

0.40 

svg 

Phase  velocity  spread  due  to  energy  spread 

1.40 

sigz 

Normalized  pulse  length 

4.2 

zO 

Normalized  Rayleigh  length 

0.42 

wO 

Normalized  mode  waist  radius 

0.65 

wmir 

Normalized  mode  radius  at  mirrors 

15 

F 

Filling  factor 

0.51 

tmir 

Normalized  mirror  separation 

19 

wbeta 

Betatron  oscillation  frequency 

0.7 

Table  5.  Dimensionless  Parameters 


The  simulation  was  used  to  mimic  the  experimental  procedure  followed  at  TJNAF 
wherein  the  electron  beam  was  tilted,  then  shifted  to  locate  the  peak  gain  for  the  given 
value  of  tilt.  For  a  given  Rayleigh  length,  the  value  of  tilt  at  which  the  peak  gain  was 
reduced  by  a  factor  of  2  was  identified.  This  process  was  repeated  for  various  Rayleigh 
lengths.  In  the  actual  laser,  the  Rayleigh  length  was  controlled  by  attaching  a  heater  to 
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one  of  the  mirrors.  The  radius  of  curvature  of  the  mirror,  which  helps  determine  the 
Rayleigh  length,  was  proportional  to  the  heater  power.  In  the  simulation,  Rayleigh  length 
is  an  input  parameter  that  is  normalized  by  dividing  it  by  the  undulator  length. 

Figure  1 1  displays  results  for  some  of  the  simulations  conducted  at  a  normalized 
Rayleigh  length  of  0.6.  It  will  be  shown  that  altering  the  Rayleigh  length  alters  the 
laser’s  sensitivity  to  shifts  and  tilts. 


Rayleigh  length  =  0.6 


“X.‘‘ 

AAAAAA*AAAa 

-0.10  0.00  0.10 
shift  (mm) 


a  tilt=0.0  mrad 
a  tilt=0.43  mrad 
a  til:=0.86  mrad 


Figure  1 1 .  Gain  versus  shift  for  a  normalized  Rayleigh  length  of  0.6 


From  Figure  11,  it  would  appear  that  the  peak  gain  for  a  given  tilt  occurs  at  or 
very  near  a  shift  of  zero.  Figure  12  shows  the  same  thing  for  a  normalized  Rayleigh 
length  of  0.4. 


Rayleigh  length  =  0.4 


~-N 

% 

V 

X 

♦  tilt  =0.0  mrad 

♦  tilt=0.43mrad 

♦  tilt=0.86  mrad 

♦  tilt=l  .28  mrad 


-0.15  -0.05  0.05  0.15 

shift  (mm) 


Gain  versus  shift  for  a  normalized  Rayleigh  length  of  0.4 
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Figure  12. 


While  the  peak  gain  for  a  given  value  of  tilt  still  appears  to  occur  very  near  zero 
shift,  Figure  12  shows  that  this  peak  is  drifting  toward  negative  shift  as  tilt  increases. 
This  effect  seems  more  pronounced  at  greater  values  of  tilt  and  smaller  values  of 
Rayleigh  length. 


Rayleigh  length  =  0.2 


shift  (mm) 


— tilt =0.0  mracf 
— *—  tilt=0.43  mrad 
— tilt=0.86  mrad 
— *—  tilt=l .28  mrad 
-Tk-  tilt=1.50  mrad 
— *—  tHt=l. 71  mrad 


Figure  13.  Gain  versus  shift  for  a  normalized  Rayleigh  length  of  0.2. 


Figure  13  is  the  same  type  of  plot,  but  for  a  normalized  Rayleigh  length  of  0.2. 
There  are  two  main  features  to  note  from  Figure  13.  First,  as  the  tilt  is  increased,  the 
peak  gain  goes  up  before  going  down.  In  Figures  11  and  12,  any  tilt  resulted  in  a 
reduction  of  gain.  In  Figure  13,  tilting  the  electron  beam  actually  increased  the  gain.  This 
effect  is  more  clearly  shown  in  Figure  14,  which  shows  tilt  versus  gain  for  a  constant 
shift  of  -.04  mm. 
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Rayleigh  length  =  0.2,  shift  =  -0.04mm 


tilt  (mrad) 

Figure  14.  Gain  versus  tilt  for  a  normalized  Rayleigh  length  of  0.2. 

Notice  that  the  peak  gain  occurs  around  a  tilt  of  1.5  mrad.  At  very  short  Rayleigh 
lengths,  tilting  the  electron  beam  can  increase  the  filling  factor,  or  overlap  between  the 
electron  and  optical  beams.  This  results  in  a  slight  increase  in  gain.  The  difference 
between  the  gain  at  zero  tilt  and  the  gain  at  a  tilt  of  1.5  mrad  is  only  2%,  so  this  is  a  small 
effect. 

The  second  interesting  feature  of  Figure  13  is  that  the  peak  gain  for  a  given  value 
of  tilt  clearly  occurs  at  a  negative  shift.  Figure  15  shows  the  relationship  between  shift 
and  tilt  with  regard  to  the  location  of  the  peak  gain.  However,  referring  back  to  Figure  13 
it  is  seen  that  the  difference  in  gain  between  the  peak  value  and  the  value  at  zero  shift  is 
less  than  1%,  so  this  is  also  a  small  effect. 
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Figure  15.  For  a  given  tilt,  this  shows  the  shift  at  which  the  peak  gain  may  be  found. 

An  additional  item  of  interest  is  seen  in  the  trend  of  peak  gain  versus  Rayleigh. 
Examining  Figures  1 1,  12  and  13,  it  is  seen  that  as  Rayleigh  length  gets  shorter,  the  peak 
gain  goes  up.  This  is  contrary  to  (4.32),  which  does  not  account  for  electron  response  to 
the  diffracting  optical  field  [27].  Since  these  diffraction  effects  are  larger  at  shorter 
Rayleigh  lengths,  it  is  not  surprising  that  (4.32)  is  less  accurate  in  this  regime. 

The  goal  of  these  experiments  and  simulations  is  to  quantify  the  relationship 
between  Rayleigh  length  and  sensitivity  to  tilt.  This  relationship  is  depicted  in  Figure  16. 
For  each  value  of  normalized  Rayleigh  length,  the  value  of  tilt  at  which  the  gain  is 
reduced  by  half  (FIWHM)  was  identified  and  plotted. 
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Figure  16.  Dependence  of  a  FEL’s  sensitivity  to  electron  beam  tilt  as  a  function  of 

Rayleigh  Length 

At  longer  Rayleigh  lengths,  the  laser  is  more  sensitive  to  tilting  of  the  electron 
beam.  At  shorter  Rayleigh  lengths,  the  optical  beam  diverges  faster,  so  the  optical  beam 
and  the  electron  beam  still  coincide,  even  when  the  electron  beam  is  tilted.  Therefore, 
short  Rayleigh  length  FELs  are  less  sensitive  to  electron  beam  tilts. 
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VII.  SIMULATIONS  TO  EXPLORE  EFFECT  OF  QUADRUPOLE 

MISALIGNMENT 


As  explained  in  section  V,  displaced  quadrupoles  alter  the  trajectory  of  the 
electron  beam.  If  a  FEL  is  to  be  used  as  a  shipboard  weapon  system,  then  it  must  operate 
in  a  vibrational  environment.  This  series  of  simulations  builds  upon  work  done  in  a 
previous  student  thesis  to  establish  vibration  control  requirements  [2],  and  was 
accomplished  in  two  stages.  First,  the  3D  FEL  oscillator  simulation  was  used  to  establish 
maximum  allowed  values  for  horizontal  and  vertical  shift  and  tilt  of  the  electron  beam. 
Then  simulations  were  conducted  using  Trace3D  and  Matlab  to  establish  peak  vibration 
amplitudes  to  ensure  the  shift  and  tilt  criteria  are  met. 

A.  TRACE3D  AND  MATLAB  SIMULATIONS 

Trace3D  was  developed  at  Los  Alamos  National  Laboratory,  and  is  an  interactive 
beam  dynamics  program  used  to  calculate  the  envelopes  of  a  bunched  beam  as  it  passes 
through  a  user  defined  transport  system  [30],  The  code  uses  transfer  matrices,  like  those 
discussed  in  Chapter  V,  to  produce  output  such  as  beam  envelopes  and  phase-space 
ellipses.  Trace3D  also  allows  the  user  to  specify  horizontal  and  vertical  displacements 
for  quadrupole  elements,  or  opt  for  random  displacements,  chosen  from  a  uniform 
distribution  with  a  user  specified  maximum. 

Due  to  the  user  interface  however,  it  is  cumbersome  to  run  numerous  simulations 
in  Trace3D.  For  this  reason,  the  transfer  matrices  were  copied  from  Trace3D  and 
incorporated  into  a  Matlab  script.  Using  loops  within  the  Matlab  script,  numerous 
simulations  were  run,  each  with  its  own  unique  set  of  randomly  generated  quadrupole 
displacements. 
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The  transport  system  segment  used  in  these  simulations  is  illustrated  in  Figure  17. 
It  is  the  same  segment  that  was  used  in  previous  simulations  of  quadrupole  misalignment 
[2],  and  is  representative  of  a  MW  class  FEL.  The  key  elements  that  require  calculations 
are  quadrupoles,  dipoles,  and  edges.  The  relevant  design  criteria  for  these  elements  are 
summarized  in  Tables  6,  7,  and  8. 


Figure  17.  Segment  of  electron  transport  system,  beginning  at  the  output  of  the 

accelerator  and  ending  at  the  start  of  the  undulator. 


Magnet 

Gradient 

(T/m) 

Location 
in  z  (mm) 

Effective  length 
(mm) 

Beam  pipe  length 
after  (mm) 

Quadrupole  1 

-6.19 

736 

150 

350 

Quadrupole  2 

6.97 

1236 

150 

350 

Quadrupole  3 

-6.16 

1736 

150 

350 

Quadrupole  4 

-1.89 

2236 

150 

175 

Quadrupole  5 

-11.7 

5491 

150 

350 

Quadrupole  6 

0.208 

5991 

150 

350 

Quadrupole  7 

-4.37 

6491 

150 

350 

Quadrupole  8 

0.0004 

6991 

150 

1706 

Quadrupole  9 

2.24 

8847 

150 

350 

Quadrupole  10 

-2.63 

9347 

150 

350 

Quadrupole  1 1 

3.85 

9847 

150 

350 

Quadrupole  12 

-7.95 

10350 

150 

300 

Table  6.  Quadrupole  design  criteria 
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The  end  result  of  introducing  displacements  for  the  quadrupoles  and  working 
through  the  transfer  matrices  is  a  set  of  four  numbers  describing  the  electron’s  position 
and  trajectory.  Horizontal  and  vertical  shift  are  specified  by  x  and  y,  while  horizontal 
and  vertical  tilt  are  represented  by  x'  and  y' .  These  numbers  could  be  entered  into  the 
FEL  simulator  to  obtain  a  value  of  extraction.  However,  characterizing  extraction  for  all 
possible  variations  of  x ,  y ,  x' ,  and  y'  would  require  tens  of  thousands  of  simulations, 
as  well  as  a  5D  plot.  Instead,  it  is  convenient  to  take  advantage  of  the  symmetry  between 
the  horizontal  and  vertical  planes. 


Magnet 

Angle  of  bend 
(degrees) 

Radius  of  curvature 
(mm) 

Location  in  z  (mm) 

Dipole  1 

-132.412 

-300 

2411 

Dipole  2 

42.412 

300 

3504 

Dipole  3 

42.412 

300 

3726 

Dipole  4 

-132.412 

-300 

4348 

Table  7.  Dipole  design  criteria 


Magnet 

Pole  face  rotation 
(degrees) 

Radius  of  curvature 
(mm) 

Location  in  z  (mm) 

Edge  1 

45 

300 

2411 

Edge  2 

-10.075 

300 

3104 

Edge  3 

48.303 

300 

3504 

Edge  4 

48.303 

300 

3938 

Edge  5 

-10.075 

300 

4348 

Edge  6 

45 

-300 

5041 

Table  8.  Edge  design  criteria 


B.  3D  FEL  SIMULATIONS 

The  3D  FEL  simulator  uses  a  linear  undulator,  which  has  a  slight  focusing  effect 
in  the  horizontal  plane.  To  achieve  symmetry,  the  undulator  pole  faces  are  curved  to 
provide  a  similar  focusing  effect  in  the  vertical  plane.  With  a  high  degree  of  symmetry,  it 
should  be  possible  to  combine  x  and  y  into  r ,  and  x'  and  y'  into  r' .  Then  extraction 
can  be  used  to  place  maximum  limits  on  r  and  r' .  These  limits  can  be  placed  into  the 
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Matlab  script  described  in  the  previous  section,  and  for  each  set  of  randomly  generated 
quadrupole  displacements  a  determination  can  be  made  as  to  whether  or  not  the  FEL’s 
performance  will  be  acceptable. 

For  these  simulations,  the  input  parameters  were  chosen  to  represent  a  MW  class 
FEL.  The  parameters  are  summarized  in  Tables  9,  10,  11,  and  12. 


Eb 

ELECTRON  BEAM  PARAMETERS 

Beam  energy  (MeV) 

97.8 

qbunch 

Bunch  charge  (nC) 

0.613 

rbx 

Beam  radius,  x  rms  (mm) 

0.1 

rby 

Beam  radius,  y  rms  (mm) 

0.1 

tb 

Pulse  duration,  FWHM  (ps) 

1.0 

prf 

Pulse  rep  frequency  (MHz) 

748.5 

lb 

Pulse  length,  FWHM  (cm) 

0.030 

gamma 

Lorentz  factor 

192 

Ipeak 

Peak  current  (A) 

613 

Iavg 

Average  current  (mA) 

459 

emitx 

Normalized  rms  x  emittance  (mm  mrad) 

2.6 

emity 

Normalized  rms  y  emittance  (mm  mrad) 

2.2 

emitl 

Longitudinal  emittance  (keV  ps) 

34 

dgog 

Beam  energy  spread  (%) 

0.08 

dthetax 

Beam  angular  spread,  x  rms  (mrad) 

0.14 

dthetay 

Beam  angular  spread,  y  rms  (mrad) 

0.11 

rho 

Beam  density  (l/cmA3) 

2.1xl014 

Pb 

Beam  average  power  (MW) 

45 

Table  9.  Electron  Beam  Parameters 
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UNDULATOR  PARAMETERS 

lambdaO 

Undulator  period  (cm) 

2.95 

N 

Number  of  periods 

20 

gap 

Undulator  gap  (cm) 

1 

Bpeak 

Undulator  peak  magnetic  field  (T) 

0.88 

Brms 

Undulator  magnetic  field,  rms  (T) 

0.62 

K 

Undulator  parameter,  rms 

1.72 

L 

Undulator  length  (cm) 

59 

Table  10.  Undulator  Parameters. 

OPTICAL  CAVITY  PARAMETERS 

S 

Cavity  length  (m) 

20 

ZO 

Rayleigh  length  (cm) 

6 

loss 

Mirror  losses  per  pass  (%) 

50 

lambda 

Optical  wavelength  (microns) 

1.6 

WO 

Mode  waist  radius,  1/e  (mm) 

0.17 

Wmir 

Mode  radius  at  mirrors,  1/e  (cm) 

2.9 

Qn 

Quality  factor 

2 

Table  11.  Optical  Cavity  Parameters. 

DIMENSIONLESS  PARAMETERS 

j 

Normalized  current  density,  linear  undulator 

83 

sx 

Normalized  beam  radius 

0.18 

sy 

Normalized  beam  radius 

0.18 

stx 

Normalized  beam  angular  spread 

0.15 

sty 

Normalized  beam  angular  spread 

0.12 

svg 

Phase  velocity  spread  due  to  energy  spread 

0.21 

sigz 

Normalized  pulse  length 

9.5 

zO 

Normalized  Rayleigh  length 

0.1 

wO 

Normalized  mode  waist  radius 

0.32 

wmir 

Normalized  mode  radius  at  mirrors 

53 

F 

Filling  factor 

0.14 

tmir 

Normalized  mirror  separation 

34 

wbeta 

Betatron  oscillation  frequency 

1.1 

Table  12.  Dimensionless  Parameters. 
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Figure  18  shows  electron  beam  shift  versus  extraction.  There  are  two  curves  in 
the  figure,  one  showing  horizontal  shift  with  zero  vertical  shift,  the  other  showing  vertical 
shift  with  zero  horizontal  shift.  Clearly,  there  is  a  high  degree  of  symmetry  between 
vertical  and  horizontal  shift. 


Figure  18.  Shift  versus  extraction  demonstrating  symmetry  between  horizontal  and 

vertical  planes. 
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Figure  19  shows  electron  beam  tilt  versus  extraction.  Again,  there  are  two  curves, 
one  showing  horizontal  tilt  and  the  other  showing  vertical  tilt,  and  again  there  is  clearly  a 
high  degree  of  symmetry  between  the  two  planes. 


0  tonsontsl  tilt  □  vertical  tilt 


Figure  19.  Tilt  versus  extraction  demonstrating  symmetry  between  horizontal  and 

vertical  planes. 

To  fully  demonstrate  the  justification  of  combining  x  and  y  into  r ,  Figure  20 
shows  a  color-coded  contour  plot  of  extraction  versus  both  horizontal  and  vertical  shift. 
The  extraction  has  been  normalized  by  dividing  through  by  the  peak  value,  so  the 
colorbar  on  right-hand  side  is  calibrated  to  overall  performance.  A  value  of  1  indicates 
that  100%  of  the  possible  output  power  is  realized.  The  radial  symmetry  of  Figure  20 
justifies  the  use  of  a  single  number  to  characterize  the  absolute  shift  off  axis. 


49 


-0. 


horizontal  shift  (mm) 


-0. 


normalized 

extraction 

n:, 
0.8 

■  -|0.6 
■  -  0.5 

■  -  0.4 

•  =  0.3 


Figure  20. 


Normalized  extraction  versus  horizontal  and  vertical  shift. 


g. 


$ 


normalized 

extraction 


I  0.6 

-  0.5 

-  0.4 

-  0.3 


2.5 

2 

1.5 

1 

0.5 

0 

-0.5 

-1 

-1.5 

-2 

-2.5 


horizontal  tilt  (mrad) 


Figure  2 1 .  Normalized  extraction  versus  horizontal  and  vertical  tilt. 

Figure  2 1  depicts  normalized  extraction  versus  horizontal  and  vertical  tilt.  Once 
again,  the  radial  symmetry  justifies  combining  x'  and  y'  into  r' .  These  combinations 
are  made  as  follows: 
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(5.24) 


r  =  ^x2  +  y2  , 
dr  dr  dx  dr  dy 

dz  dx  dz  dy  dz  ’ 

r'  =  i(x2  +  v2)  1/2  (2x)x'  + x 2  +  /) '/_  (2v)/  , 


(5.25) 


In  general,  it  is  reasonable  to  consider  anything  more  than  a  20%  degradation  of 
performance  unacceptable.  In  this  case,  a  20%  degradation  of  the  peak  extraction,  which 
is  2.77%,  with  an  average  electron  beam  energy  of  45MW,  results  in  an  output  optical 
power  of  <  1  MW.  Therefore  it  is  reasonable  to  choose  cutoff  values  of  r  and  r'  based 
on  a  normalized  extraction  of  0.8.  The  un-normalized  data  for  Figures  20  and  21  is 
presented  in  Tables  13  and  14.  Performance  of  80%  of  the  full  potential  corresponds  to 
an  extraction  of  -2.2%.  The  data  in  Tables  13  and  14  has  been  color-coded  with  red 
indicating  acceptable  performance  and  blue  indicating  unacceptable  performance.  From 
these  tables  we  see  it  is  reasonable  and  conservative  to  require  r  <  0.4  mm  and 
r'  <  1.1  mrad. 


Table  13.  Extraction  versus  horizontal  and  vertical  shift. 
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Table  14.  Extraction  versus  horizontal  and  vertical  tilt. 


Figure  22  shows  the  percent  of  the  time  the  above  criteria  will  be  met  versus  the 
maximum  allowed  value  of  displacement,  assuming  a  uniform  distribution.  Figure  23 
shows  the  percentage  versus  the  standard  deviation  of  displacement,  assuming  a  Gaussian 
distribution. 


Uniform  distribution  of  quadrupole  displacements 


Figure  22.  Reliability  versus  peak  displacement  for  a  unifonn  distribution. 
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Gaussian  distribution  of  quadrupole  displacements 


Figure  23.  Reliability  versus  standard  deviation  of  displacement  for  a  Gaussian 

distribution. 

From  Figures  22  and  23  we  can  see  that  to  ensure  an  acceptable  level  of  optical 
output  power  100%  of  the  time,  uniformly  distributed  quadrupole  displacements  must  be 
kept  less  than  ~  0.02  mm,  or  20  microns.  For  a  Gaussian  distribution  of  quadrupole 
displacements,  100%  reliable  requires  the  standard  deviation  to  be  kept  less  than  .007 
mm,  or  7  microns. 
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VIII.  CONCLUSION 


It  has  been  shown  that  in  order  to  function  as  a  weapon  system,  a  FEL  must  have 
an  output  optical  power  on  the  order  of  1  MW.  The  perfonnance  can  be  reduced  due  to 
shipboard  vibrations.  It  is  has  not  been  determined  what  the  exact  output  power  of  the 
eventual  weapon  will  be,  nor  is  it  known  beyond  a  rough  estimate  how  much  power  is 
required  to  successfully  engage  an  incoming  missile.  Therefore,  the  acceptable  amount 
of  performance  degradation  can  only  be  estimated.  Using  conservative  estimates,  it  has 
been  found  that  quadrupole  alignment  on  a  vibrating  ship  must  be  kept  to  within 
tolerances  on  the  order  of  10  microns.  This  is  consistent  with  previous  estimates,  and 
within  the  capability  of  active  alignment  systems  already  in  use  [31,  32].  Such  a  system 
will  be  required  to  employ  a  FEL  onboard  a  ship,  but  these  tolerances  do  not  comprise  an 
insurmountable  engineering  problem. 
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